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We combine field observations, microcosm, stoichiometry, and molecular and stable isotope techniques to 
quantify N2O generation processes in an intensively managed low carbon calcareous fluvo-aquic soU. All the 
evidence points to ammonia oxidation and linked nitrifier denitrification (ND) being the major processes 
generating N2O. When NH4"^-based fertilizers are applied the soil will produce high N2O peaks which are 
inhibited almost completely by adding nitrification inhibitors. During ammonia oxidation with high NH4"^ 
concentrations (>80 mg N kg"') the soil matrix wiU actively consume oxygen and accumulate high 
concentrations of N02~, leading to suboxic conditions inducing ND. Calculated N2O isotopomer data show 
that nitrification and ND accounted for 35-53% and 44-58% of total N2O emissions, respectively. We 
propose that slowing down nitrification and avoiding high ammonium concentrations in the soil matrix are 
important measures to reduce N2O emissions per unit of NH4"'" -based N input from this type of intensively 
managed soil globally. 

Nitrogen-fertilized intensively managed agricultural soils represent the largest global anthropogenic 
source of the potent greenhouse gas N20''^. Understanding the N2O production processes is the pre- 
requisite for developing targeted mitigation strategies. Most previous studies have been conducted under 
laboratory conditions and it is still a challenge to connect these results to real field situations'''''. 

N2O emissions from soils are produced by a diversity of microorganisms^ through at least three main micro- 
bially-mediated processes nitrifier [nitrification (NN), nitrifier denitrification (ND) and denitrifier denitrifica- 
tion (DD)] although abiotic reactions are also known to be an important N2O source under certain 
circumstances". The occurrence of each N2O production process and its contribution to total N2O emissions 
depend on prevailing soil conditions such as soU oxygen conditions (water status), microbiaUy available C and N, 
and pH, which are governed by edaphic and climatic conditions and agricultural management practices'"''. 
Although many approaches have been employed for source partitioning of soil N2O emissions under laboratory 
conditions, none of the current methods has convincingly quantified the N2O sources under real field situa- 
tions'''". However, a combination of existing methods provides great potential to narrow the uncertainties 
associated with N2O source partitioning in the field. 

Our laboratory and field experiments conducted over a number of years investigating intensively managed low 
carbon calcareous soils support the assertion that nitrifying prokaryotes are responsible for a large share of the 
emitted N2O" '^ in marked contrast to many European studies on agricultural soils with higher readily decom- 
posable carbon and moisture contents and with relatively low pH values in which N2O is produced mainly by 
denitrification. However, we still have no unequivocal evidence that the high peaks of N2O emission induced by 
applying NH4*- or urea-based fertilizer to the soils studied are derived from the first steps of nitrification (NH4* 
to NH2OH to NO2") or from successive steps of ND (NO2" to NO to N2O to N2). Indeed, one of our earlier 
studies strongly suggested an important role of ND in N2O production". We therefore hypothesize that the 
transition from aerobic (NN) to suboxic (ND) conditions might be largely responsible for the N2O peaks after the 
application of NH4*- or urea-based fertilizer. 

Here, we report a series of pot incubation experiments to study the response of N2O emission to soil water 
content, NH4^ concentrations, and varying N and C supply. Moreover, a robotized incubation/analytical system"" 
was employed to study the stoichiometric relationship between N rate, gaseous dynamics (N2O, NO, N2, O2 and 
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CO2), and the response of these gases to the type of N and C sources. 
Two field experiments were conducted to investigate the isotopic 
composition of N2O and the composition and abundance of the 
functional genes of nitrification (amoA) and of nitrifier denitrifica- 
tion or denitrification (nirK and nirS). The experiments were carried 
out in different years and under a range of conditions using the same 
soil type at two sites on the North China Plain in order to show that 
the phenomenon and results are congruent and support each other. 

Results 

N2O emission vs. soil water content, C, and N supply. Total cumu- 
lative N2O emissions were significandy (P < 0.05) related to soil water 
filled pore space (WFPS) in urea or urea plus straw treatments, 
with maximum emissions at around 70% WFPS. In the treatment 
without urea, total cumulative N2O emissions were very low, and only 
slightly affected by soil moisture, even at high soil NOs^-N status 
(31.8 mg N kg"') (Pot Incubation Experiment 1, PIEl) (Fig. la). 
These observations are very different from previous reports based 
on European agricultural soils with high indigenous soil NOs^-N 
status". We speculated that the large amount of N2O in 70% WFPS 
was produced from the nitrification process in urea or urea plus straw 
treatments because the course of the N2O emission was coincident 
with urea hydrolysis and the nitrification process''. Measurements of 
soil NH4*-N and NOs^-N status at the end of the incubation period 
partly explain the conclusion (fi^om substrate bases as side evidence) 
that the occurrence of high N2O emissions at 70% WFPS could be 
attributed mainly to nitrification and ND, or coupled nitrification- 
denitrification (CND) (see also Discussion section). The NH4* 
contents at higher moisture levels (90 and 110% WFPS) were 
significantly higher than at lower moisture levels (30, 50 and 70% 
WFPS), whereas their NO3" contents were significandy lower 
(Table S5, P < 0.05) irrespective of whether nitrogen and straw 
were added to the soil. This indicates that nitrifiers might dominate 
N-transformation processes when the soil water content is below 70% 
WFPS but denitrification might dominate N-transformation pro- 
cesses when the soil moisture exceeds 90% WFPS. 

Cumulative N2O emissions exhibited a significant positive linear 
correlation with NH4+-N input rate (R^ = 0.9898, P < 0.001) under 
70% WFPS (Pot Incubation Experiment 2, PIE2) (Fig. lb). NH4+-N 
input rate almost completely explained the variation in total N2O 
emissions and the emission factor (here, the slope of the linear func- 
tion) was 0.57% in this incubation system, very similar to previous 
field experiments in the summer maize season'-*''^. Soil NH4*-N, 



NOs^-N andN02~-N contents at days 1, 3 and 18 during the incuba- 
tion period showed that NH4*-N was rapidly (within 1-3 days) 
nitrified to N02"-N and then NOs^-N and that high N02"-N accu- 
mulated when NH4+-N > 80 mgNkg"' (Table S6). The nitrification 
process corresponds well with the time course of N2O emission (Fig. 
S2). The NO3" content at the end of incubation was closely linearly 
related to NH4+ input rate (R' = 0.9922, P < 0.001) (Fig. S3). This 
further confirms that N2O emissions were dominated by NH4"'' 
transformation process. 

Application of a nitrification inhibitor (NI), either DMPP or DCD, 
reduced the cumulative N2O emissions by 99.2% and 97.1%, respect- 
ively, compared to urea alone. They almost completely blocked N2O 
emission in urea or urea plus straw treatments (Pot Incubation 
Experiment 3, PIE3) (Fig. Ic). Similar to PIE2, soil NH4+-N was 
rapidly nitrified to NOs^-N in urea and urea plus straw treatments 
and this process slowed down when using NI (Table S7). 
Furthermore, NI application precluded N02"-N accumulation so 
that no more N2O was produced (Table S7; Fig. S4). 

Stoichiometry of N transformation and gas fluxes. We performed 
an experiment with different urea-N concentrations using a robotized 
incubation/analytical system" to investigate the relationship between 
the emissions of various gases and O2 consumption (Robotized Incu- 
bation Experiment 1, RIEl). Cumulative N2O emission increased 
exponentially (P < 0.001; Fig. 2a) and CO2 production increased 
linearly (P < 0.001; Fig. 2b) with increasing concentration of urea- 
N. The O2 consumption slope of the linear relationship was only 1.45 
which lies between the theoretical slopes of pure nitrifier nitrification 
and pure nitrifier denitrification (Fig. 2c), suggesting that the N2O 
production originated from a combination of these processes. The 
significant increase in CO2 may be attributed to the hydrolysis of urea 
and the concomitant release of CO2. The dynamics of O2 consump- 
tion, CO2 production, and N2O emissions were synchronized with 
each other and this demonstrates well the close coupling between 
urea hydrolysis, nitrification, and N2O production. 

N2O emission patterns were significantly different between high 
concentrations of urea-N (320 mg N kg"') and low concentrations 
(<80 mg N kg"') which likely reflected the different N2O produc- 
tion processes involved. NO and N2 emissions were generally low 
within the concentration range 0-80 mg N kg"' but gradually 
increased at the end of incubation in the highest urea treatment 
(320 mg N kg"'). This can be attributed to the depletion of O2 by 
nitrifiers and the induction of denitrification (Fig. S5). At the end of 
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Figure 1 | Cumulative N2O emissions in response to soil conditions in pot microcosm incubation experiments, (a), N2O emission as a function of soil 
water fUled pore space (WFPS). CK denotes control treatment. U denotes applied urea at a rate of 40 mg N kg"'. US denotes application of urea 
with comminuted wheat straw, (b), N2O emission as a function of NH4*-based fertilizer (ammonium sulfate) input rate at 70% WFPS. (c). Effect of 
nitrogen (urea), carbon source (comminuted wheat straw), and nitrification inhibitors (DMPP and DCD) on soil N2O emission at 70% WFPS. 
DMPP and DCD denote 3,4-dimethylpyrazole phosphate and dicyandiamide, respectively. Error bars in plots (a-c) are standard deviation of the mean 
(n = 3 replicates). 
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Figure 2 | Stoichiometric relationships between N2O, CO2, O2 and urea-N input rate, (a), N2O emission as a function of urea-N input rate, (b), CO2 
emission as a function of urea-N input rate, (c), O2 consumption as a function of urea-N input rate. In plot (c), NN and ND denote nitrifier 
nitrification and nitrifier denitrification, respectively. Note that O2 consumption in the highest urea input treatment (U320) was not included due to the 
development of anaerobic conditions in the later stages of incubation and the units of x- and y-axis in plot (c) are different from plots (a & b). 



the incubation period the O2 concentration in the 320 mg N kg"' 
treatment was 2.36%, low enough to trigger denitrification'". Our 
automated incubation experiments demonstrate that denitrification 
can proceed immediately once the oxygen consumption has lowered 
O2 to a low level. The importance of O2 conditions on the denit- 
rification process outweighs the easily decomposed carbon limitation 
reported in a previous study'^. 

To further explore the kinetics of N gas production under condi- 
tions favoring nitrification or denitrification, we performed another 
robotized incubation/analytical experiment {RIE2) with different 
conditions in terms of NH4*-N, NOs^-N, organic carbon amend- 
ment, NI amendment, O2 concentration, or combinations of these 
factors (Fig. 3). N2O concentration increased to 3.9 ppm at the 
beginning of the incubation with NH4*-N (Fig. 3b) and was blocked 
by simultaneous application of a nitrification inhibitor (Fig. 3c). 
However, there was virtually no N2O emission when NOs^-N 
was applied alone or together with glucose (Fig. 3d,e). In contrast. 



both the N2O concentration (from 61.1 to 1307.5 ppm) and the 
NO concentration (from 16.0 to 1834.9 ppm) increased sharply 
and showed clear peaks at the initial stage of incubation when 
NOs^-N was applied together wdth glucose under anaerobic condi- 
tions (Fig. 3f). NO and N2 emissions were not detected in the other 
treatments except when N03"-N was applied together with glucose 
under anaerobic conditions, in which NO, N2O, and N2 production 
patterns typical of denitrification occurred. Under aerobic condi- 
tions, NO, N2O, and N2 emissions were not elevated even when 
glucose was applied with NO3" -based fertilizer and this supports 
our conclusion from RIEl that oxygen is the dominant factor inhib- 
iting denitrification rather than carbon limitation. We therefore con- 
clude that 70% WFPS was too low to trigger denitrification in our 
soil. However, if we withdraw O2 from the system the typical denit- 
rification process will proceed. At 70% WFPS aerobic nitrification 
dominated soil N transformations and N2O was the major gaseous 
product. 
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Figure 3 | Kinetics of gaseous N2O, NO, N2, O2, and CO2 in robotized analytical/incubation systems, (a), Control treatment (CK). (b). Ammonium 
sulfate (AS), (c), AS with nitrification inhibitor (AS -I- NI) 3,4-dimethylpyrazole phosphate (DMPP). (d). Calcium nitrate (CN). (e), CN with 
glucose (CN -I- C). (f), without oxygen from the system of treatment CN -I- C (CN -I- C-O2). Soils were incubated at 70% WFPS. Error bars in plots (a-f) 
are standard deviation of the mean (n = 4 replicates) . Note that the right and left y axis scales have large differences. In some cases the error bars are smaller 
than the symbols. 
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Isotopic fingerprints of soil-derived N2O. N2O generated by nitrifi- 
cation is more depleted in '^N and '"O in relation to precursor N 
compounds compared with N2O derived from denitrification". 
Furthermore, the distinction in '""N site preference (SP:difference 
in §'^N between the central and peripheral N positions of the 
asymmetric N2O molecule) between nitrification and denitrifi- 
cation differentiates these two N2O sources^". We established a 
field experiment to study the N isotopic characteristics of soil N2O 
derived from three treatments, namely fertilized with (NH4)2S04 
(AS), Ca(N03)2 (CN) or Ca(N03)2 plus compaction of the soil, 
glucose addition, and increased irrigation (creating soil environ- 
mental conditions favoring denitrification, CN -I- C). The N2O 
emission patterns were detailed in a previous paper" and are 
summarized in Fig. S6. The characteristics of 8"*0-N20, §'^N-N20 
and A8'^N-N20 are detailed in Fig. 4a,b,c. 

In the AS, 5"*0-N20 (referenced to standard mean ocean water, 
SMOW) increased with time after fertilization (27.9 to 36.4%o). The 
time gradient of S^'O was 2.2%o d"'. In the CN, S'^O of N2O 
increased during the first three days after fertilization (36.6 to 
45.0%o) with a time gradient of 4.2%o d"' and levelled off on the fifth 
day with a 8"*0 of 44.2%o. In the CN -f- C, §'"0 of N2O increased 
from 26.6 to 43.4%o during the first three days after fertilization (the 
data on the fifth day were lost) with a time gradient of 8.4%o d"' 
(Fig. 4a). The rapid enrichment of '"O may be attributed to the 
intensive consumption of NO3" and the subsequent enrichment of 
'"O in the remainder of the NO3" in the soil as shown by the high N 
fertilizer N2O loss. In this treatment N2O emission accounted for 
52.9% of the input N during the 5 days following fertilization 
(Fig. 4f)^'. Although the overlaps of S'^O from the soil favoring 
denitrification with the other two treatments and the complicated 
O sources of N2O derived from nitrification and denitrification make 
it difficult to partition N2O sources based solely on S'^O^^'^^, our 5"*0 
for N2O in the AS was within the upper end of the reported S'^O for 



N2O produced in aerobic environments (13 to 35%o)^'', suggesting 
that nitrification was the predominant N2O source. 

The S'''N (referenced to standard atmospheric N2) of N2O 
increased with time after fertilization in all three treatments. In the 
AS, S'^N increased from - 33.5 to - 7.5%o. The time gradient of 5'"N 
was 6.5%o d"' and the calculated A6'^N of N2O on the first sampling 
day was -58.7%o (= difference between 5'=N of N2O and 5'=N of AS 
of which the 5'=N was 25.2%o). In the CN (8'=N = 9.0%o), 5'^N 
increased from —7.0 to 1.6%o with a time gradient of 2.2%o d~' 
and the calculated A8'^N on the first sampling day was — 16.0%o. 
In the CN -I- C, 5'^N increased from —1.3 to 15.2%o with a time 
gradient of 8.3%o d"' and the calculated AS'^N on the first sampling 
day was — 10.3%o (Fig. 4b). The S'^N of the third sampling event was 
6.2%o higher than that CN and we attribute this to the enrichment of 
'^N in the remainder of the soil NO3" by the same mechanism of 
enrichment as discussed for '"O above. 

Our AS'^N values on the first sampling day in the AS agree very 
well with the reported AS'^N values for nitrification ( — 111 to 
— 35%o)^^"^'', indicating that nitrification was the main source of 
N2O. The A5'''N values in the CN are within the range of reported 
AS'^N values for denitrification, and in the CN + C they were close to 
the lower end of the range (—45 to — 13%o)^^'^'', indicating that denit- 
rification might be an important source of N2O in the former. We 
attribute the relatively low enrichment factor in the latter to the 
reduction of N2O resulting in 5'''N enrichment of the remaining 
N2O. In the CN + C treatment a plot of S'^O versus S'^N'"'"' of 
N2O gave a slope of 0.89 revealing that N2O reduction might have 
occurred to some extent (Fig. 4e)^^. 

In the AS, SP remained constant with a mean value of 15.4%o 
(ranging from 12.7 to 17.1%o). The CN SP showed an increasing 
trend and the mean value was 15.0%o (ranging from 11.4 to 
23.5%o). In the CN + C, SP decreased slightly over time with a mean 
value of-2.4%0 (ranging from -5.4 to0.1%o) (Fig. 4c).OurSP ofthe 
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Figure 4 | Isotopic fingerprints and emissions of NjO. (a), S'^O. (b), b'^N''"^. (c), 5'=N site preference (SP). (d), 8"'0 (vs. SMOW) versus SP. (e), 5"*0 
versus S'^N''""'. (f), cumulative N2O emissions. AS, CN, CN + C, and CK denote amendment with ammonium sulfate, calcium nitrate, and calcium 
nitrate -I- glucose -I- compaction of soil, and control, respectively. Error bars in plots (a), (b), (c), and (f) are standard deviation of the mean (n = 4). Note 
that data points are individual observations of the treatments in plots (d) and (e). 
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soil favoring denitrification (CN + C) was very consistent with the 
SP of N2O produced by bacterial denitrification in pure culture stud- 
ies (0 to -lOXoY"'^". In the AS, SP was lower than typical results 
reported for the NH20H-to-N20 pathway from pure cultures of 
nitrifiers (33%o) but was higher than the SP of denitrifiers or nitrifier 
denitrification (0 to — 10%o)^"'^''. Similar results were observed in our 
earlier soil incubation study in which the same soil was incubated at 
45% WFPS and amended with AS and the observed SP was 18.4%o". 
Therefore, both NH2OH oxidization and NO2" reduction contrib- 
uted substantially to the release of N2O from the AS. Furthermore, 
the positive correlation between S"*0 and with a slope of 0.29 

(P < 0.00 1 ) of the N2O also indicates that nitrifier denitrification was 
a significant N2O source (Fig. 4e)^'^'^°. In addition, the lack of a rela- 
tionship between S'^O and SP indicates a constant relative contri- 
bution of NH2OH oxidization and NO2" reduction to the gross N2O 
production in AS and CN + C (Fig. 4d) because the SP values for 
these two N2O production processes are constant while their §'"0 
values are variable due to the different O sources and to O exchange 
with H2O as well as the variable 5"*0 of the O sources"''^. 

In the present study the SP value of soil-derived N2O in the AS 
treatment indicates that both NH2OH oxidization (NN) and NO2" 
reduction (including ND and DD) processes contribute substantially 
to total N2O production. Moreover, our molecular analysis also 
shows that it was nitrifiers that dominated NH4* transformation 
(see paragraph below). Evidently, the additional N2O derived from 
NO2" reduction in the AS treatment compared to the CN treatment 
can be attributed to the ND process. Therefore, we calculate the 
fertilizer-induced N2O generating processes in the AS treatment 
according to the following equations: 

Fas = FAS(Totai) — Fck = Fas(nn) + Fas(nd) + Fas(dd) ( 1 ) 

FaS(DD) ~FcN(Total) — FcK (2) 

FaS(NN) = FAS(Total-NN) — FcK(NN) (3) 

FcK(NN)~FcN(NN) (4) 

FAS(Total-NN) = FAS(Total) X ^NN' FcN(NN) = FcN X ^NN (5) 

fNN = (SPbulk — SPdemt)/(SPnit " SPdenit) (6) 

where F is cumulative N2O emissions, AS is ammonium sulfate, CN 
is calcium nitrate, Total is the gross N2O including both fertilizer 
derived and background emissions, CK is the control treatment, NN 
is nitrifier nitrification, ND is nitrifier denitrification, DD is denitri- 
fier denitrification, and Total-NN is NN produced N2O emissions 
including both fertilizer and background derived. f^N is the fraction 
of N2O derived from the NN process, SPbuik is the measured SP value 
of N2O, SPnit and SPjemt are the SP values for nitrification and 
denitrification processes, respectively. There are two assumptions 
needed to proceed to the source partitioning. First, the DD derived 
N2O was equivalent to Ca(N03)2 fertilizer derived N2O in the treat- 
ment receiving CN only. This assumption is reasonable because most 
of the applied NH4* was transformed into NO3" as shown in our 
microcosm studies, and denitrification thus contributed to the total 
N2O emission of the AS treatment. Second, that NN process derived 
N2O in the control was equivalent to that of NN in the treatment 
receiving Ca(N03)2 only. Then we adopt flux- weighted SP value and 
the above mass-balance equation with two end-members to calculate 
the relative contribution of NN, ND, and DD to the gross (NH4)2S04 
and Ca(N03)2 derived N2O of the fertilizer treatments, respectively. 
In the two end-members equation we adopt 33%o as the SP for 
nitrification and 0%o or — 10%o for denitrification in two rounds of 
calculation"'^"'^". Finally, we show that NN, ND, and DD account for 



35-53%, 44-58%, and 2-9% of (NH4)2S04 derived N2O emissions, 
respectively. 

N20-generating microorganisms and functional genes. To further 
understand the linkage of functional genes associated with N2O 
production we investigated the community structure and the 
abundance of the ammonia oxidizing bacteria (AOB) and 
denitrifiers during the crucial periods of N2O emission. The 
treatments were zero-N (Control, CK), NH4*-based N (AS), 
NH4+-based N + DMPP (AS + NI), NOs^-based N (CN), and 
N03"-based N -I- glucose -I- more water (CN + C). The N2O emis- 
sion patterns in all treatments were very similar to those previously 
reported''-'^- (Fig. S7). 

A high abundance of amoA gene copy numbers for AOB was 
found in AS treatment but DMPP clearly suppressed amoA gene 
copy numbers and also impacted nirK/nirS gene copy numbers to 
some extent (Fig. S8a). CN + C treatment increased nirK/nirS gene 
copy numbers considerably but effects from the CN treatment were 
negligible (Fig. S8b,c). Community structures of AOB in AS and AS 
+ NI treatments in the pivotal period of N2O emission were similar 
to that in the CK treatment at each sampling time (Fig. S9a). There 
were no changes in the nirK and nirS gene community structure 
under various N-fertUizer treatments during all crucial periods of 
N2O emission (Fig. S9b,c). 

The phylogenetic tree revealed that all AOB sequences were 
affiliated with Nitrosospira or Nitrosospira-like species (Fig. SIO). 
The results indicate that Nitrosospira species are ubiquitous in the 
soil studied and this is consistent with the results of other studies 
conducted in this region'^ and the findings in most soO media world- 
wide^'''. The dendrogram of deduced amino acids shows five major 
clusters of nirK sequences (I to V) (Fig. Sll). The great majority of 
clones were closely related to nirK of Rhizobiales as well as clones 
which belonged to cluster II""^. However, these n/rK-harboring deni- 
trifiers do not belong to Rhizobiales because nirK phylogeny is 
incompatible with 16S rDNA phylogeny""". The nirS tree was divided 
into four clusters (clusters I-IV) (Fig. S12). The major nirS clones 
were grouped in clusters I and IV and were primarily related to nirS 
of Rhodobacteraceae and Rhodocyclaceae^^ . Only two nirS (102 bp) 
clones were grouped in cluster II; while the remaining clones 
belonged to cluster III closely related to the nirS of Micrococcaceae 
and Burkholderiaceae^". 

Conceptual model of N2O generation and regulation. Fig. 5 pre- 
sents a summary conceptual model of N2O generation and regulation 
in this calcareous soil. Ammonia oxidation is the first N2O -generating 
process. It is regulated by NH4^-N concentration and soil moisture 
and temperature and can be efficiently inhibited using nitrification 
inhibitors. Ammonia oxidation occurs in aerobic conditions but 
consumes oxygen significantly and accumulates NO2" temporarily, 
leading to suboxic conditions that induce a secondary nitrifier- 
denitrification process to generate N2O with high NH4* concentra- 
tions in the soil matrix. Denitrifier denitrification is less active under 
normal field conditions with high O2 concentrations in the light 
textured soil which has low moisture during most of the year and 
low available carbon. Due to the very low denitrification activity the 
soil profile normally accumulates large amounts of N03"-N which 
presents great leaching potential to the shallow groundwater when 
surplus N is present in intensively managed high N-input agricultural 
systems"*'"'''''"'. We conclude that N2O emissions are due mainly to the 
activities of nitrifiers in this calcareous fluvo-aquic soil with 35-53 
and 44-58% of N2O produced by nitrification and nitrifier 
denitrification, respectively. 

Discussion 

Collectively, all of our evidence points to NN and ND as the dominant 
N2O generation processes in this intensively managed calcareous 
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fluvo-aquic soil with DD playing a minor role. Our results are further 
supported by a recent report showing that ammonia oxidation path- 
ways (NN and ND) are important sources of N2O under low oxygen 
availability following urea or ammonium sulfate application to loam, 
clay and sandy loam soUs (that also contain ample nitrate) by using 
advanced '^N-'^O isotope and NH3 oxidation inhibition approaches'". 
The weak contribution of DD to N2O emission was substantiated by 
the weak response of N2O emissions to soil NO3" content in our pot 
incubation experiments and in the robotized incubation experiments 
and the very low N2O emission factors of the soil receiving NO3" 
based fertilizer in the field" '*. Although DD can proceed in anaerobic 
microsites in the soil matrix according to numerous reports''"''''\ the 
weak contribution of DD to N2O emission in the soil studied was 
ultimately attributable to the lack of anaerobic microsites for trigger- 
ing the reduction of NO3" and thus the subsequent sequential reduc- 
tion of N oxides by denitrifiers and the associated N2O production. 
The weak contribution of DD to N2O emission on the North China 
Plain has also been reported in a previous study and has been attrib- 
uted mainly to the low content and availability of soil organic carbon 
for denitrifiers'^. However, our robotized incubation experiments 
demonstrate that DD can proceed immediately once the oxygen con- 
sumption has lowered O2 to a low level. The importance of the soil 
redox conditions on the process of DD outweighs that of dissolved 
organic carbon (DOC) limitation as demonstrated in RIE2 within 
which N2O emission responded weakly to the simultaneous addition 
of NO3" and glucose under the normal oxygen levels. Therefore, the 
foremost influence of DOC on DD in this region might be its indirect 
effects on creating anaerobic soil conditions by consimiption of soil 
O2 and the formation of soil aggregates within which anaerobic 
microsites tend to resulf*^'**. 

The dominant contribution of N2O production pathways, includ- 
ing NN and ND, to N2O emissions from NH4* or urea-based fert- 
ilizer was substantiated by the observed intensive and proportional 
response of N2O emissions to soil NH4* and NO2" contents in our 
laboratory and field studies. The overwhelming importance of the 
NH4* oxidation pathway on N2O emissions in the present soil has 
also been verified by two of our earlier '^N-labeled laboratory incuba- 
tion investigations in one of which NH4^ oxidation accounted for 
88% of total N2O emissions at 45% WFPS", and in the other NH4+ 



oxidation accounted for 85 to 95% and 83 to 94% of the total N2O 
emissions at 40% and 60% WFPS'^, respectively. These two studies 
also demonstrate that the relative contribution of nitrification and 
denitrification to gross N2O emissions was quite constant across this 
relatively large range of water content in these soils. Surprisingly, the 
relative contributions of nitrification and denitrification reported in 
these two studies agree well with our present field study in which 
NH4* oxidization and DD accounted for nearly 90 and 10% of total 
(NH4)2S04-derived N2O emissions, respectively. However, neither 
of our former '^N-labeled studies discriminated between the contri- 
butions of NN and ND to N2O derived from NH4* oxidization. 

It has long been debated how microbial community structure or 
abundance and ecosystem functions are coupled. Our study suggests 
that environmental factors (mainly soil N content) affected the 
abundance of AOB and denitrifiers harboring nirK and nirS genes 
positively and linkages occur between N2O emissions and the abund- 
ance of active AOB and denitrifiers harboring nirK, and nirS genes 
rather than their community composition. There was a striking vari- 
ation in the abundance of AOB and denitrifiers in our different 
treatments and the shifts coincided with N2O emission patterns, 
suggesting that shifts in active abundance play a relatively important 
role in determining N2O emissions under different N, C, and water 
management strategies. 

The potential for nitrifiers to proceed with denitrification and 
associated N2O production has been recognized for a long time'"'''''"' 
and was recently verified to be a prevalent trait of the beta-proteo- 
bacterial AOB such as Nitrosomonas europaea and Nitrosospira spp. 
in pure culture investigations*' and to be a significant N2O source in 
soil incubation investigations'' '"''*''. In our field study all of the AOB 
sequences were affiliated with Nitrosospira or Nitrosospira-like spe- 
cies, allowing us to elucidate the processes of N2O production as well 
as the associated SP values for NN, ND, and DD. 

We did not identify the response of ammonia-oxidizing archaea 
(AOA) to N fertilization and it has been reported that AOA are 
ubiquitous in soils but they do not respond to NH4* oxidization 
and N2O production in intensively managed agricultural soils^'. 
The contribution of archaea to the SP value of the associated N2O 
is not yet clear^''™'^'. A recently reported SP value of N2O from 
marine AOA is virtually identical to that of AOB"**. Up to now the 
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soils where AO A have a significant impact on NH4* oxidization 
seem to be mainly confined to acid types^^. One study conducted 
in our region also revealed that Nitrosospira-like AOB were dom- 
inant over AOA in oxidizing NH4*^'. 

We did not have the direct evidence to show that the N02^ pro- 
duced by ammonia oxidization will not be utilized by denitrifiers in 
CND process, but is utilized instead by nitrifiers in ND. However, it 
produces different isotopic and isotopomeric values, which gave rise 
to the separation in the first place for isotopic studies. In addition, our 
soil molecular results showed that only the nitrifiers respond to the 
NH4*-based fertilizer application and the denitrifiers only respond 
to the NO3" -based fertilizer together with glucose application. These 
results reinforce our conclusion that the reduction of NO2" in our 
soil should mainly through ND pathway while the contribution of 
DD pathway might be minor. Some previous studies have concluded 
that CND is not a separate process^ but instead comprises rapid 
denitrification of the nitrification byproduct and end product 
(NO2" & NOa")^''. This process has been found in paddy soUs^'', 
wastewater^* and soil aggregates'^ when favorable environmental 
conditions for both nitrification and denitrification prevail, but the 
soU-climatic conditions in our study are very favorable for nitrifica- 
tion processes". 

The conventional wisdom, based on research on soils with higher 
readily decomposable carbon, higher moisture contents and rela- 
tively low pH values in Europe, is that denitrification by denitrifiers 
is the dominant soil N transformation process leading to N2O emis- 
sions. In contrast, the soil type in our studies (calcareous fluvo-aquic 
soil) is widely distributed among agricultural soils globally and is 
frequently associated with semi-arid or sub-humid continental cli- 
mates, and typically has high pH, low available C, aerobic and low 
moisture conditions. The dominant effect of nitrifiers and the 
importance of their denitrification function and the consequentially 
large amounts of N2O emissions may be phenomena that are wide- 
spread globally, at least in similar edaphic-climatic conditions. Our 
understanding of N2O production processes can contribute to the 
development of efficient mitigation strategies for these land areas, i.e. 
to slow down nitrification and avoid high ammonium concentra- 
tions in the soil matrix. 

Methods 

Study sites. Two sites were used in the present study, one located at Dongbeiwang 
town, Beijing {40^^00' N, 116 12' E) and the other at Quzhou county, Hebei Province 
(36''52' N, 115" 10' E), China. Fields at the two sites have typical intensively managed 
agricultural soils on the North China Plain where a winter wheat -summer maize 
rotation is the dominant crop production system and urea and NH4^-based fertilizer 
are the most commonly applied nitrogen fertilizers. The soils are calcareous Eluvo- 
aquic soils (calcareous Cambisols according to the FAO Classification). SoU 
properties in the upper 20 cm at the Dongbeiwang and Quzhou sites were: texture 27/ 
57/16 and 39/58/3 (sand/ silt/clay, %); bulk density 1.37 and 1.34 g cm'^ CaCOs 
content 52.5 and 58.2 g kg"'; DOC 51 and 43 mg C kg"'; pH (1 : 2.5, soil/water) 7.9 
and 8.2; EC 0.1 8 and 0.16 mScm"'; total N 1.2 and 0.7 g kg"'; organic matter 2 1.3 and 
14.7 g kg"', respectively. 

Pot incubation experiments (PIE). Three PIEs were conducted to investigate the 
response of N2O emission to soil water content {PIEl ), to NH4"^ input rate {PIE2), and 
to varying nitrogen and carbon source supply {PIE3). Surface soils (0-20 cm) used for 
PIEs were collected from Quzhou site at the beginning of April, 2010. InPIEi,it was a 
completely factorial design with three factors, namely control, urea (applied at a rate 
of 40 mg N kg"' dry soil), urea plus straw (applied at a rate of 1.0 g C kg"' dry soil), 
and within each factor five classes of soil water content, 30, 50, 70, 90, 1 10% WFPS 
were included (Table SI). Incubation was conducted for 15 days. lnPIE2, nine classes 
of N input rate, 0, 10, 20, 40, 80, 160, 240, 320, 400 mg N kg"' dry soil were selected 
and (NH4)2S04 was the selected N fertilizer. Incubation was conducted at a water 
content of 70% WFPS for 18 days. In PIE3, soils were amended with and without N 
fertilizer (urea, applied at a rate of 40 mg N kg"' dry soil) and wheat straw (applied at 
a rate of 1.0 g C kg"' dry soil) and nitrification inhibitors (DMPP and DCD) (Table 
S2). Incubation was conducted at a water content of 70% WFPS for 16 days. All 
incubations were conducted aerobically and in the dark at 20^C with three replicates 
of each treatment. N2O emissions were measured periodically by gas chromatography 
(GC). Detailed information about the operation of PIEs and N2O measurements can 
be found in the SI. 



Robotized incubation/analytical experiment (RIE). Detailed information about the 
robotized incubation/analytical system has been given by Molstad et aV^. In brief, 
the robotized system comprises an incubation system and a gas analysis system. The 
headspace air of the incubation flasks was replaced with a mixture of helium (He) and 
oxygen (O2) at the desired proportions. Gas in the headspace of each flask was 
sampled every 8 hours by CTC GC-PAL (delivered by Agilent) and analyzed by gas 
chromatography (Agilent GC 7890A) for the concentrations of N2O, N2, 02, and CO2 
and by NO analyzer for NO. Two RIEs were conducted to investigate the 
stoichiometric relationship between N input rate and gaseous (N2O, NO, N2, 02 and 
CO2) flux (RIEl) and the response of these gas fluxes to the type of N and C sources 
and inhibitors {RIE2). Soils used for RIEs were collected from the Quzhou site in AprU 
201 1 . In RJEl , six classes of N input rate, 0, 20, 40, 80, 160, and 320 mg N kg" ' dry soil 
were adopted and urea was selected as N fertilizer. In RIE2, six treatments were set up, 
including control (CK), amended with (NH4)2S04 at a rate of 40 mg N kg"' dry soil 
(AS), amended with nitrification inhibitor DMPP on the basis of AS (AS + NI), 
amended with Ca(N03)2 at a rate of 40 mg N kg"' dry soU (CN), amended with 
glucose (1.0 gCkg"' dry soil) on the basis of CN (CN + C), and anaerobic incubation 
on the basis of CN + C (CN + C-O2). RIEs were conducted in a thermostatic water 
bath at 20"C and at a soil water content of 70% WFPS with four replicates for 256 and 
272 h for RIEl and RIE2, respectively. 

Isotopic fingerprint experiment. The experiment was conducted in the field at 
Dongbeiwang site in the maize season in 2007, based on Field study 2 reported in Ju 
et al}^ of which three treatments (AS, CN, and CN + C) were selected to measure the 
N2O isotope signatures of soil gases (Table S3). In summary, treatment AS was a 
typical N fertilizer management practice for maize on the North China Plain, i.e. 
application of NH4"^-based N fertilizer ((NH4)2S04 in the present study, 150 kg N 
ha"', the same rates as other fertilizer treatments) followed by irrigation. Treatment 
CN was amended with Ca(N03)2 and received the same amount of irrigation as 
treatment AS. Treatment CN + C was set up to create conditions favoring 
denitrification, i.e. applying N03"-based fertilizer (Ca(N03)2 in the present study) 
and glucose, compacting the soil, and increasing irrigation. Each treatment was 
replicated four times. The methods for collecting and storing gas samples for isotope 
analysis are introduced in detail in the SI. Isotopic signatures of N2O (S'^N*"^, S'^N", 
5'^N'^, S'^O) were determined by isotope ratio mass spectrometry (IRMS, Thermo- 
Finnigan, Bremen, Germany) and the results are expressed in 5-notation as per mil 
(%o) deviations from the isotopic ratios of atmospheric N2 or standard mean ocean 
water (SMOW). For 6'^N" and S'^N**, pure N2O was used as reference gas which was 
calibrated in the laboratory of N. Yoshida at the Tokyo Institute of Technology. The 
difference between the '^N isotopomer ratios of the central and peripheral N of N2O 
(S'^N"* - S'^N"*) is referred to as '^N site preference (SP). Detailed information on 
isotopic analysis of N2O and calculations were given by Well et a/.". 

Microorganisms and functional gene analysis. Soils were collected from a field 
experiment conducted at Quzhou site in the maize season in 2009'*. The field 
experiment consisted of five treatments (Table S4), namely (1) control treatment 
without N input (CK), (2) amended with NH4"^-based fertilizer (AS), (3) amended 
with nitrification inhibitor DMPP and NH4+ -based fertilizer (AS + NI), (4) amended 
with N03"-based fertilizer (CN), (5) amended with glucose and N03"-based fertilizer 
(CN + C). (NH4)2S04 and Ca(N03)2 were used as the NH4"'- and N03"-based N 
fertilizers (250 kg N ha"'), respectively. The N fertilizers were split into two 
applications: 100 kgNha"'atthethird-leafstage(June27'^) andl50 kgNha"'atthe 
tenth-leaf stage (July 29""). We directly irrigated 20 mm after each fertilization. Soil 
was sampled on two periods of strong N2O emissions (July T' and August 3'"'^) and two 
periods of low N2O emissions (August 7* and September 26*). A DNA-based 
approach with real-time PCR and T-RFLP analysis was used to assess the abundance 
and community structures of the AOB and denitrifiers and focused on the amoA gene 
for AOB and the nirK/nirS gene for denitrifiers. Detailed information about microbial 
and functional gene analysis and N2O emission measurements can be found in SI. 

Statistical analysis. The primary data were treated using Microsoft Excel 
spreadsheets. The data shown in this paper were calculated by averaging three or four 
replicates. The cumulative N2O, soil mineral N, gene copies and T-RFs among the 
different treatments were tested by analysis of variance and mean values were 
compared using SAS statistical software (Version 8.2; SAS Institute, Inc., Gary, NC) to 
calculate least significant difference (LSD) at the 5% level. Total N2O fluxes and WFPS 
and NH4"^-N application rate were analyzed by regression using Sigmaplot 11.0 
software (Systat Software Inc., Erkrath, Germany). 
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